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1 Introduction 

One of the most interesting questions of the classical differential geometry which has ap- 
peared at studying of semi-Hamiltonian systems of hydrodynamical type is the description 
of the surfaces admitting not trivial deformations with preservation of principal direc- 
tions and principal curvatures. Then the number of essential parameters on which such 
deformations depend, is actually equal to number various local Hamiltonian structures 
of corresponding system of hydrodynamical type [6]. Such local Hamiltonian structures 
are determined by differential-geometrical Poisson brackets of the first order (see [10]). 
In the same paper a bi-Hamiltonian structure of the system of averaged one-phase so- 
lutions of KdV have been considered. Later multi-Hamiltonian structures of systems of 
hydrodynamical type were studied in [12], [13]. The given work is devoted to the de- 
scription of bi-Hamiltonian structures in language of orthogonal curvilinear coordinate 
nets. The alternative approach (by the method of an inverse scattering transform) has 
been used in [16]. In [5] it has been shown, that the description of pairs compatible 
differential-geometric Poisson brackets of the first order is equivalent to the solution of 
the over-determined system on rotation coefficients j3 ik of orthogonal curvilinear coordi- 
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nate nets 

diPik + dkPki + J2m^i,k PmiPmk = 0) * 7^ &j (1) 

where ^(r 1 ) are some functions, <9fc = d/dr k , k = 1, 2, iV. The general case r/j(r l ) 7^ 
const and particular case ^(r 1 ) = q = const were described in [5], where it has been 
shown, that (1) is an integrable system. In that specific case ^(r*) = q = const the 
system of equations 

^iAfc + <9fc/^fci + 12m^i,k PmiPmk = 0) * 7^ &j (2) 

when iV = 3 was reduced to three pairwise commuting non-evolution equations (see [5]) 
& t = chqV 1 - 9x> ?j/t = shg^l - gg, ^ = -^(1 - q*)(l - q*), 

which are nothing but two modified Sin-Gordon equations and degenerated twice modified 
Sin-Gordon equation, respectively. In this paper generalization of this system (i.e. %(r l ) ^ 
const) is presented for bi-Hamiltonian structure where one metric is flat, and the second 
one has curvature of co- dimension 1. 



2 Flat case 

Compatibility condition of the linear system in partial derivatives 

diH k = f3 lk H h i^k (3) 
is the same as compatibility condition for conjugated linear system 

9i^ k = PkiiPr, i k - ( 4 ) 

This yields the nonlinear multi-dimensional integrable iV-wave system (see for instance, 
[4], [9]) 

diP jk = PjiPik, i^J^k. (5) 

Any pair of solutions of the first linear system H^, generates hydrodynamic type 
system integrable by the generalized hodograph method [14] 

rr(2) 

r\ = ' ri, • = 1, 2, N, 



written in Riemann invariants r 1 . Particular solutions tp^ of the second linear system (4) 
determine densities a k and fluxes c k of conservation laws 

a\ = d x c k (&), 

where 

N N 

da k = ^t ] H$dr m , dc k = ^g^Wdr™ 

m=l m=l 

Existence of the linear differential operator of the first order connecting solutions of both 
linear systems (see for instance, [13]) 

H i = dilpi + J^^mi^m (6) 

is equivalent to existence of local Hamiltonian structure 

a] = d x [g ik C 



da k 

where flat coordinates a k are determined by "null" solutions of the second linear system 

= d^\ k) + J2P mi ^\ fc = l,2, ... ,N, 

and non-degenerate symmetric metric 

TV 

ds 2 = ( H mdr m ) 2 = g jk da j da k , 

m=l 

which is constant in flat coordinates a k 

N 



9 3k = Yj^S^S = const • 



m=l 



Existence of such operator (6) imposes the restriction on rotation coefficients of orthogonal 
curvilinear coordinate nets (3 ik 

9if3 ik + d k f3 ki + Y PmiPmk = 0, i^k, (7) 

well-known in differential geometry as the Gauss equation, in above case it fixes the metric 
of zero curvature. 

It is easy to see that iV-wave system is invariant under the transformation R l = R % {r l ), 
Vi 1/2 (r l ), where ^(r l ) and R % (r l ) are arbitrary functions of a single variable. If 
Riemann invariants are fixed (R l (r l ) = r l ), then "flatness" condition (7) is no longer valid, 
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— 1/2' 

then solutions of the conjugate problem (4) transform in accordance with ip i = ip i fi i [r % 
if fJ>i(r l ) = R 1 ' (r l ), then a metric preserves "flatness" 



N N 
, 2 



m=l m=l 

and solutions of the conjugate problem (4) are the same 

TV N 



m=l m=l 

Thus, infinite number of metrics 

N 



ds 2 = Y.^ rm ) ( H $dr m y 



m=l 



are connected to each system of hydrodynamic type and only for finite number of values 
/ij(r*) (no more than N + 1, see [7]) the "flatness" condition is satisfied. In this paper 
the case is considered, when just two distinct values fi^r 1 ) determine flat metrics. Since 
coefficients of the diagonal metric Hf are determined up to multiplication on an arbitrary 
function of corresponding Riemann invariant, then the first flat metric always can be cho- 
sen so that one of values ^i{r l ) can be fixed to 1. Compatible pairs of local Hamiltonian 
structures determined by such metrics studied for instance in [16]. Such pairs are com- 
patible if their arbitrary linear combination also is local Hamiltonian structure. Thus, the 
integrable nonlinear system (see [5]) 

diP jk = PjiPik, i^j^k, (9) 

(A + RWPik + (A + R k )dkhi + l$ik + hi) + E ( A + Rm )P™Pmk = 0, i*k. 
is a compatibility condition of the linear system 



= (A + + ) + Emd X + Rm )Pmi^l\ k = 1, 2, ... , N, 



(10) 



where 3 = d/dR\ j3 ik (H) = (3 ik (r)/R l '(r % ), = i^j'K A is an arbitrary constant of the 
flat metric g n (r) = (A + R l )g %% . 

Thus, we have proved following 

Lemma: The reduction of the iV-wave system (5) 

diPik + dkPki + ^2mjti,k PmiPmk = 0; * 7^ k, 

Vl (rW lk + Vk(r k )8 k p kt + Ir/iWfa + \i k {r k )P lk + £ m ^. fc V m (r m )P mi P mk = 0, i ? k, 
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where rj^r 1 ) is a second solution /^(r 1 ), is resulted in (9) by scaling H^r) = rf^r^H^R), 
d/dr 7, = rjKr^d/dK 1 , r/^r 1 ) = fj^R 1 ) = R\ if r/^r 1 ) ^ C{ = const. 

Thus, in general case bi-Hamiltonian structures of hydrodynamic type systems are 
described by solutions of the integrable system 

r f 9^ tt + r fc «9 fc /^ + §/3 ifc + \(3 ki + £ m ^. fc r™/3 mi /3 mfe = 0, i^k. 
Its spectral problem (10) 

= (A + r*)d t 4 k) + l4 k) + E m ^(A + k = 1, 2, iV 

is nothing but iV commuting systems of ordinary differential equations with respect to 
every independent variable (Riemann invariant r l ). 

Remark: This multi-dimensional spectral problem has a mult i- dimensional analog of 
Wronskian: N(N + l)/2 "first integrals" are constraints 

N 

sT = E( A + rm )V&tyS? = ^nst, (11) 

m=l 

which are nothing but metric coefficients in flat coordinates of mixed local Hamiltonian 
structure. 

Following G. Darboux for further consideration we should determine "first integrals" 
of the Gauss equation (see (7)). 

Definition: A scalar potential of rotation coefficients of conjugate curvilinear coor- 
dinate nets V is a such function, whose mixed second derivatives are 

Vik = PikPki, i + k - 

Definition: A vector potential of rotation coefficients of conjugate curvilinear coordinate 
nets S k is a such function, whose mixed first derivatives are 

diS k = (3 ik d k (3 ki , i ^ k. 

G. Darboux have proved (see [3]), that the Gauss equation, written via rotation coefficients 
of conjugate curvilinear coordinate nets (7), has N "first integrals" 

Si + \Y,0 2 rm-<r% (12) 

where rii{r l ) are arbitrary functions ("integration factors"). Similar zero curvature con- 
dition 

VitfWu. + Vk(r k )d k f3 kt + ^)f3 lk + \i k {r k )(3 lk + £ vJr m )(3 mi (3 mk = 0, i^k 

mj^i,k 
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for metric Tj i {r' l ')g n has N "first integrals' 1 



where ki(r l ) are another arbitrary functions. 



2.1 Particular case, rj^r 1 ) = C{ = const 

This particular case was considered in [5]. The spectral problem 

= (A + C % )d^f + Zmd X + C m)^t\ k = 1, 2, 

determines the nonlinear system (2) 

It easy to see that the consequences of (12) and (13) 



CiS'i + 2 X^m^i C m(3 m i ~ &i 



for particular case ^(r 1 ) = c, = const are the constraints 

5^(cm - Ci)f3 2 mi = 2[ki(r l ) - CiUiir 1 )}. 

Thus, if iV = 3, then above system 

(c 2 - ct)(3 2 21 + (ca - Ci)p 2 31 = hir 1 ), 

( Cl -c 2 )p 2 12 + (c 3 -c 2 )pl 2 = l 2 (r 2 ), 

( Cl -c 3 )(3 2 13 + (c 2 -c 3 )f3 2 23 = Z 3 (r 3 ), 

can be parameterized 



/3 2 i = \ coshu, (3 31 = \ sinhii 





(r 1 ) 


c 2 


- Cl 


k 


(r 2 ) 


Cl 


- c 2 


\ 


(r 3) 


Cl 


- c 3 



/?i2 = V coshu, P 32 = \ sinhu, 

y c i — c 2 y c 2 — c 3 

fl / ^(r 3 ) / / 3 (r 3 ) . 

P13 = A/ coshw;, (3 23 =\l- -smhw. 



ill 


r 1 ) 


Cl 


- c 3 


y 




c 2 


- c 3 




(r 3) 


c 3 


- c 2 



Then the 3-wave system has the form 



-diw = t -. 2 ° 3 -coshw, — d^u = a/7 ~T7~ rsinhu>, 



\/h(r l ) V (ci-c 2 )(ci-c 3 ) v^sP) V (ci -c 3 )(c 2 - c 3 ) 



o / C 2 - C 3 1 / d - C 3 

=OiW = W7 77 r-smhw, — o 2 u = 4/7 77 rsinnu, 



\A( rl ) V (ci - c 2 )(ci - c 3 ) v// 2 (r 2 ) V (c 2 - Ci)(c 2 - c 3 ) 

<9 2 u> = a /t 1 3 rcosht;, — 83V =4/7 1 2 rcoshw. 



^(r 2 ) V (ci - c 2 )(c 3 - c 2 ) v// 3 (r 3 ) V ( c i ~ c s)(c 2 - c 3 ) 

Introducing new independent variables (scaling Riemann invariants) 



P = V (* - c 2 )'(c?- o) J • /W ~ )dTl ' 



ci - c 2 



(ci - c 3 )(c 2 - c 3 ) 
finally one can obtain the system (see [5]) 



d p w = coshw, d s u = sinhw, 
d p v = sinh u, d q u = sinh v, 
dgiu = coshu, d s v = coshw. 



For instance, eliminating field variables v and w (v = arcsinh-Ug, w = arcsinhu s ), one 
can obtain two modified sinh-Gordon equations 



u pq = sinhuwu 2 + 1, u sp = coshwyw 2 + 1 



and the degenerated twice modified sinh-Gordon equation ([1]) 

U sq = ^(u 2 + l)( U 2 + l). 

Such equations can be obtained as well as with respect to other field v and w 



v P g = smhvJvl + 1, v sq = cosIiva/i; 2 - 1, v sp = J (v 2 p + \){v 2 s - I) 



w 



sq 



coshwJw 2 -!, w sp = sinhw Jw 2 - 1, w pq = J (w 2 - l)(w 2 - 1) 



The substitution 

z — v + w — w + ln[w q + \Jw 2 — 1] 
transforms these equations into 



z sq = sinh z, z ssp = z s J zl + z 2 - z p . (14) 



The first one is well-known Bonnet equation (the sinh-Gordon equation), the second one 
is the next commuting flow from a hierarchy of the potential modified KdV equation (see 
[15]) 

1 



3 

*T %SSS 2^ S ' 



Thus, a spectral problem for a triple of the modified sinh-Gordon equations has a form 





where ip 1 = y/c 2 — c^ip, ip 2 = V c 3 ~ c i^5 ^3 — V c i ~ a spectral parameter is 

(A + ci)(c 2 - c 3 ) 



(A + c 3 )(ci - c 2 ) 



Remark: Each above linear system can be reduced to well-known spectral problem 2x2 
for the modified sinh-Gordon equation (see [1]) by transformations described in [17]. 
This becomes obvious if take into account that all three above linear systems 3x3 have 
constraint (see (11)) 

(ijj 2 + (1 — C)^ 2 ~~ 4> — const . 

It means that the functions ip, x/j, tp are "squares of eigenfunctions" for corresponding 
spectral transform 2x2 (see [17]). Also, every above linear system 3x3 can be written 
as a scalar spectral problem of the third order. For instance, the third linear system is 
equivalent to the third order linear equation 

^Psss """^ss — (C + sinh 2 w + w^)'i/j s + [(( + sinh 2 w)—^- — 3w s sinh w cosh w]ip = 0. (15) 

w s w s 

At the same time, the well-known fact is (see for instance, [8]) that the Yajima-Oikawa 
hierarchy (the Yajima-Oikawa system also is known as the long-short resonance) is asso- 
ciated with the spectral problem 

14 = (16) 

where 

L = d 2 s — a + a\d~ 1 a 2 . (17) 

Substituting (17) in (16) and comparing with (15), one can obtain expressions for the 
coefficients 

a = sinh 2 w + w 2 , a\ = w s , a 2 = w ss — sinh w cosh w. 



S 



That means, that the sinh-Gordon equation is embedded not only to the KdV hierarchy 
but also to the Yajima-Oikawa hierarchy. In other words: the pseudo-differential Manin- 
Sato operator 

L = d s + A d; 1 + A 1 d; 2 + ... 

associated with the KP hierarchy has finite-component reductions [9] 

JV 

L n = a: + A , n d:~ 2 + ... + A n . %n + 



fc=l 

where N is an arbitrary natural number. Thus the KP reduction determined by operator 

L = d 2 — (sinh 2 w + w 2 ) + w s d~ 1 (w ss — sinh w cosh w), 

describes three-component bi-Hamiltonian structures of hydrodynamic type systems. 

Remark: One can introduce field variable c = z p , then re-scale independent variables 
d s — > ed s , d p — > e _1 9p, then a second equation from (14) has a form 

2a„_a c + £2 ° ss 



e d s c = d. 



^/c 2 + e 2 c 2 s 



A limit of this equation with respect to the parameter e when e — > yields 

d c = d \—- —} 

which is nothing but again the sinh-Gordon equation 

b sp = sinh b, 

where ln c = b. 

Remark: The substitutions 

z 1 = w ± v, z 2 = V ± u, z 3 = u±w 

connect solutions of the sinh-Gordon equation 

"I 1 O O Q Q 

z sq = sinh z , z pq = sinh z , z sp = cosh z . 

2.2 General case, rj^r 1 ) = r l 

In general case rj^r 1 ) = r % the spectral problem 

= (A + r*)dM k) + l 2 4 k) + EmM X + rm )P™^t\ k = 1, 2, ... , N, 
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determines the nonlinear system 

di(3 jk = Pjifiik, i^j^k, 

diftik = r i- r k [— \{fiik + Pki) + S m ^j,fc( r ~ rm )PmiP m k\i i k. 

It is easy to see that consequences of (12) and (13) yield another "first integrals' 1 
For instance, if iV = 3, then above constraints 



(r 1 - r 2 )/^ + (r 1 - r 3 )/? 2 



31 



(r 2 - r 1 )/^ + (r 2 - r 3 )/? 2 2 



^3- 



describe a bi-Hamiltonian structure determined by the flat metrics g %% and r % g %% . 
Using parametrization 



An 
/5 32 







13 



^ 1 ry 2 



coshw, /3 31 



Vi 



^3 ^* 1 



ry 2 ryi> 



cosh f, /3 



12 



^ 1 ^ 2 



V 3 



,y 3 ry 



- cosh iw, 



28 



V* 



ry 2 ry^> 



sinh u, 



sinh t>, 



sinh w, 



the system (5) transforms into 



d 3 u 



13 



2(r 3 — r 1 



sinh u sinh iw 



2(r 3 -r 1 ) 



^•3 1 



7- cosh iu cosh u 



cosh w sinh u, 



ry 2 ^«3 



^ 1 ,y«2^ ^,y«3 ry 1 ^ 



^ 1 ry 2 



^^3 ^" 1 ^ ^ryQ 7""^^ 



Vi cosh u, 



Vs sinh to, 



<9i-u 



r, 



12 



1/ ^2 



2 ( r l _ r 2) 



2/ * 1 



2(t-1 _ r 2^ 



cosh u cosh v + 



sinh v sinh w + 



^ 2 ^*3 



^ ^ 1 ^*2 ^ ^i^«3 ^ 1 ^ 



— v Vi sinhu, 



^*3 ry 1 



^i^il ry 2 ^ ^ ^ 2 3 ^ 



Vo cosht;, 



sinh f cosh tt, 



^* 1 ^* 2 
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d 2 w = —7-7^5 rr cosh v cosh to + W— — — — vV^sinht;, 

2(r^ — r d J y (r 1 — r z ){r z — r 6 ) 



VW/V2 . , . , , / r 1 - r 2 



^ = -7T7^ 7- sinh w sinh t; + W— — — VV^coshw, 

^ ( y 7™ ) 1/ ( T* 7" ) ( y* y ] 



V23 = -5 5- cosh v sinh it), 

where Vik = PikPki- Each above triple of equation depends from third independent variable 
as a parameter, which can be eliminated by shift along both other independent variables. 
Thus, every above system can be written in form 

d x w = ^7 x ^ v sinh u sinh w + , I ^ -a/KcosIi-u, 
2(x -y) V x ( x ~ V) 

d v u = ^ V ^ \ cosh io cosh m + , / -"^4 sinh to, 

\/VxVy 

V xv = — cosh to sinh u. 

x-y 

3 Co-dimension 1 

In previous section a flat case is considered, where relationship between two conjugate 
linear problems (3), (4) is determined by the linear differential operator of the first order 
(6). In general case of arbitrary curvature matric (8) such relationship becomes nonlocal 
(see for details [11]) 

M M 

H % = + (3 m ^ m + ^5> n iff (18) 

m^i fc=ln=l 

where c^a^ = , Ekn are constant symmetric non-degenerate matrix, and a number 

of particular solutions of (3) is determined by the Gauss equation (see [11]) 

M M 

mj^i,j k=l n=l 

Thus, in general case the description of pairs of nonlocal differential operators (18) (similar 
to the flat case) is equivalent to the system 

difijk = PjiPiki » ^ 3 + k, 

M M 

8$^ + djPji + J2m&,j PmiPmj = Efc=lE,=l^^ H j , * ^ 3, 
. . MM _ 

r% diPij + ^^jPji + 2@ij 2@ji ^2m^i,j ^PmiPmj = Y2k=l Y2n=l £ knH i Hj , % ^ j. 
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In case of non-flat metric t] i (r z )g u the Gauss equation 

1 1 M M 

<l,(r')<), i irllj{ rJ)<hi ir ^{,(r : )!, r ^(y)i, r ^ Vm (r m )P mi f3 mj = J^ 6 ^ ^ 

m^i,j k=l n=l 

has N "first integrals" 

1 1 .J M M 

+ ^ + 2 E ^ = 2EE^^ fc) ^ n) + k ^> 

mj^i fc=ln=l 

where /^(r*) are some functions. 

In this section we restrict our consideration on two simplest cases: a flat metric g w 
and a metric g n of co-dimension 1, when g l% = Cig %% or g u = r l g n . In the first case the 
system 

+ djPji + Em^j AnAtf =0, % ± j, (19) 
c idiPij + c jdjPji + Em^ij C mPmiPmj = HiHj, % ^ j, 

is a result of a compatibility condition of the linear system (see [16]) 

diipj = fljiipi, i ^ j, 

= (A + C^dilpi + Em^i( A + C m)P mi 1p m + Hid, 

did = Hupi. 

In the second case the system 

diPij + Ojpji + Em^ij PmiPmj =0, 1 7^ j, (20) 

is a result of a compatibility condition of the linear system (see [16]) 

diipj = fljiipi, i ^ j, 
= (A + r*)^ + |Vi + E m ^(A + r m )/3 m ^ m + H,a, 
did = Hupi. 
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3.1 Particular case, rj^r 1 ) = q = const 
As well as in the flat case the nonlinear system (19) 



has the constraints 



CiSi + | ^2im+i C mPmi = \Hf + k i{ r% )- 



Thus, one has 



£(* - c m )PL + H? = 2[cMr l ) ~ Hr% 
When A = 3, then above system has a form 

(c 1 - C2 )^ 1 + ( Cl -C 3 )/32 1+ /f2 = 

(c 2 - Ci )P 2 12 + (c 2 -c 3 )P 2 32 + H^ = Z 2 (r 2 ), 
(c 3 - ci)/3? 3 + (c 3 - c 2 )^ 3 + Hi = h(r 3 ). 

Introducing new independent variables (re-scaling Riemann invariants) 
and new dependent functions 



Hi — y/hRi, H 2 = yhS2, H 3 = V/3P; 



3, 



Ci - c 3 V Ci - c 3 V Ci - c 2 

23 = ~i\ — - — Pi, P12 = H — - — S 3 , (3 32 = \ — - — Si, 
V c 2 - c 3 V ci - c 2 V c 2 - c 3 
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the system (5) can be written in form pairwise commuting hyperbolic systems 
dpPr = iR 3 P 2 , d s R l = -JL=P 3 R 2: 



VI - A 2 " 



d p P 2 = -iP 1 R 3 -^—^P 3 R 1 , d s R 2 = iR 3 P x - J^— R 1 P 3 , 



A 



d p P 3 — R\P 2 , d s R 3 — —iP\R 2 , 



d p S 1 


— R2S3, 




d q Ri = 


AS 2 R 3 , 


d p S 2 


— -^RiS 3 , 




d q R 2 = 


SiR 3 , 


d p S 3 


= —R 2 S\ - 


- -^R\S 2 , 


d q R 3 = 


~- —S\R 2 — AS 2 R±, 


d q Pi 


= tS 3 P 2 + 


iVl - A 2 P 3 S 2 , 


d s S 1 = - 


-s 3 p 2 y/T ^s i ^ 


d q P 2 


= -iS 3 Pt, 




d s S 2 = 





d q p 3 = -iVT=&s 2 Pi, d s s 3 = p 2 s 1 . 

where 

a 2 c 2 -c 3 



ci - c 3 



p( + p 2 i + Pi = 1, R( + R z 2 + Ri = i, s( + si + si = i. 



(21) 



Each above system is nothing but a some reduction of the famous Cherednik model of 
chiral fields related with hierarchy of the Heisenberg magnet (see [2]). 

3.2 General case, rj^r 1 ) = r l 
In general case the nonlinear system (20) 

9iP ik = ^[H.H, - l((3 ik + (3 kl ) + Em^ k ( rk - rm )PmiPmkl * ^ k. 

has a set of "first integrals" 

Si + \ J2m^i @ mi = n i{ r% )i 

r*S t + \Vi + \ Em?i r m /3L = \m + Hr 1 )- 

Then we have 

Vi + ]>> m - r^L = H? + 2fa(r i ) - amir*)]. 
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Thus, if TV = 3, then the system of 18 equations (21), written in form of three pairwise 
commuting hyperbolic flows 



&1023 = 021013, 9 3021 = 02303H 

d\H 3 = @ 13 Hi, 9 3 ifi = /3 31 iJ 3 , 

an — H-iHq, _ gi3+^31 _ r 2 -r 3 a a p, o _ _HjJh _i_ Pw+Psi _ r x -r 2 a a 

u ll J 13 ~ r i_ r 3 2(r 1 -r 3 ) r i_ r 3 / u 21^23' U 3H31 ~ r i_ r 3 "1" 2(r 1 -r 3 ) r i_ r 3 A^21/ u, 23' 

<9l/? 32 — 031012) ^2/^31 — 031021-, 

d 1 H 2 = Pi 2 Hi, d 2 Hi = 21 H 2 , 

ft R — H 1 H 2 _ ffl2±ffgl I r 2 -r 3 g n n n _ H 1 H 2 _i_ ^12+^21 _ ^-r 3 /3 /3 

u ll J \2 ~ r i_ r 2 2(r 1 -r 2 ) ^ r 1 -r 2 -^31^32' ^2/^21 ~~ r l_ r 2 T 2 (r 1 -r 2 ) r i_ r 2M3l/- , 32) 

^2/^13 — '120 '23' ^3/^12 — ^13^32' 

•92-^3 = 023^2, d 3 H 2 = 32 H 3 , 

p> ft _ H^Hz _ _ r'L-r 3 a a ft R H2H3 1 023+032 1 ■r 1 -r 2 a a 

U 2H23 ~ r 2_ r 3 2(r 2 -r 3 ) r 2_ r 3 M12M13' ty 3/ J 32 ~~ r 2_ r 3 ~r 2(r 2 -r 3 ) ~ 1 ~ r 2 -r 3 ^12/ c 'l3' 

has three "first integrals" 

(r 1 - r a )/& + (r 1 - r 8 )/& + = Vi, 

(r 2 - r 1 )^ + (r 3 - r 8 )^, + flf = V 2 , 

(r 3 - r x )0\ 3 + (r 3 - r 2 )0i 3 + R\ = V 3 . 



Since this example and for higher dimensions and co-dimensions presence of "the first 
integrals" any more does not lead to an essential reduction of number of the equations 
in each of such hyperbolic systems. For example, at N = 3 from 18 equations in case of 
two flat metrics where ratios of corresponding diagonal coefficients are various constants, 
- there are 6 equations of the first order; if ratios of corresponding diagonal coefficients 
are some functions of Riemann invariants, - there are 9 equations of the first order; if one 
of metrics is flat, and the second one is of co-dimension 1, and ratios of corresponding 
diagonal coefficients are various constants, - there are 12 equations of the first order; and 
at last, if ratios of corresponding diagonal coefficients are some functions, - there are 15 
equations of the first order (see explanation below), - in all other more complicated cases 
(any dimension of space N and any co-dimension M) the number of equations in system 
in involution will be more the of use of "first integrals", than in initial system of the 
equations. For example, in last above-stated example - the system has 18 equations. Use 
of 3 constraints on a first step decreases number of the equations up to 12. However, 
as constraints contain the first derivative of scalar potential V of rotation coefficients of 
conjugate curvilinear coordinate nets, it is necessary to add 3 more equations expressing 
second derivative from V through rotation coefficients of conjugate curvilinear coordinate 
nets. Thus, actually already in this example it is visible, that above system contains the 
same number of equations (as any equation of the second order is equivalent to a pair of 
the first order equations). 



4 Conclusion 

In this paper it has been shown, that the problem of classification of semi-Hamiltonian sys- 
tems of hydrodynamical type with such additional and important properties as a flatness 
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of the diagonal metrics or metric's curvature of co-dimension 1 is reduced to the solu- 
tion of such well-known integrable systems (by inverse scattering transform) as a negative 
flow of the Landau-Lifshitz equation (i.e. the Cherednik model) and generalization of the 
Ernst equation well known in Gravity and General Relativity (here a spectral problem 
depends explicitly of independent variable). 
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